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ABSTRACT 

We present new multi-wavelength observations of the dwarf Seyfert 1 galaxy POX 52 in order to 
investigate the properties of the host galaxy and the active nucleus, and to examine the mass of its 
black hole, previously estimated to be ~ 10^ Mq. Hubble Space Telescope ACS/HRC images show that 
the host galaxy has a dwarf elliptical morphology (A// — —18.4 mag, Sersic index n = 4.3) with no 
detected disk component or spiral structure, confirming previous results from ground-based imaging. 
X-ray observations from both Chandra and XMM-Newton show strong (factor of 2) variability over 
timescales as short as 500 s, as well as a dramatic decrease in the absorbing column density over a 9 
month period. We attribute this change to a partial covering absorber, with a 94% covering fraction 
and Nh = 58l^;^ x lO^i cm , that moved out of the line of sight in between the XMM-Newton 
and Chandra observations. Combining these data with observations from the VLA, Spitzer, and 
archival data from 2MASS and CALEX, we examine the spectral energy distribution (SED) of the 
active nucleus. Its shape is broadly similar to typical radio-quiet quasar SEDs, despite the very low 
bolometric luminosity of Lboi = 1-3 x 10'*'^ ergs s~^. Finally, we compare black hole mass estimators 
including methods based on X-ray variability, and optical scaling relations using the broad H/3 line 
width and AGN continuum luminosity, finding a range of black hole mass from all methods to be 
Mbh = (2.2 - 4.2) X 10^ Mq, with an Eddington ratio of Lboi/^Edd « 0.2 - 0.5. 
Subject headings: galaxies: individual (POX 52) — galaxies: active — galaxies: dwarf — galaxies: 
nuclei — galaxies: Seyfert — X-rays: galaxies 



I. INTRODUCTION 

Stellar and gas dynamical studies have proven to be the 
most secure method of measuring the masses of super- 
massive black holes in nearby galaxies, and black holes 
with masses of ~ 10^ - 10^ Mq have now been de- 
tected in a few dozen galaxies (see Ferrarese & Ford 2005 
for a recent review). With just a few exceptions, how- 
ever, stellar-dynamical and gas-dynamical techniques 
for black hole mass measurement cannot be applied to 
Seyfert galaxies and quasars, because most luminous ac- 
tive galaxies are too distant for their black hole's grav- 
itational sphere of influence to be resolved. Most black 
hole mass estimates in broad-lined active galactic nuclei 
(AGNs) are based on indirect methods that rely on scal- 
ing relations resulting from rev erberation-mapping data 
(jKaspi et al.ll2000l : lBentz et al. 2006). Using these meth- 
ods, it has recently become possible to search for active 
galaxies having extremely low-mass black holes, in order 
to extend our understanding of black hole demographics 

^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. These 
observations are associated with program # 10239. 

^ Based on observations obtained with XMM-Newton, an ESA 
science mission with instruments and contributions directly funded 
by ESA Member States and NASA. 

^ Department of Physics & Astronomy, 4129 Fredrick Rcincs 
Hall, University of California, Irvine, CA 92619-4575; thorn- 
toc@uci.edu, barth@uci.edu 

* The Observatories of the Carnegie Institution of Washington, 
813 Santa Barbara Street, Pasadena, CA 91101 

Department of Physics, McGill University, 3600 rue University, 
Montreal, QC, H3A 2T8, Canada 

® Department of Astrophysical Sciences, Princeton University, 
Princeton, NJ 08544 

^ Hubble Fellow and Princeton-Carnegie Fellow 



to objects having black hole masses of Mbh ^ 10^ Mq. 

The best nearby example of a black hole in the < 10^ 
Mq mass ran ge is that in the dwarf Seyfert 1 galaxy 
NGC 4395 (.FiliDpenko fc SargentJl98"llFiliDDenko fc Hoi 
120031 ). NGC 4395 is an essentially bulgeless, late- 
type dwarf galaxy with type 1 Seyfert ch aracteris- 
tics such as broad emission lines (jFilippenko fc Sargen^ 
[1989; Filippcnko et al. 1993 !) and ra pid X-ray variabil- 
ity (Iwasawaet al..,2000; .Shih et all (20031: iMoran et all 
Located at a distance of only ~ 4.3 Mpc 
(|Thim et al.l[2fl0l . it is an excellent target for multi- 
wavelength observations despite the tiny luminosity of 
its central engine (Lboi ~ 6 x lO'*'^ ergs s"'^; Moran et al. 
1999). A recent ultraviolet (UV) reverberation-mapping 
measurement found M bh = (3.6 ± 1 .1 ) x 1 0^ Mq 
(jPeterson et all |2005[) . [G reene IFlfol (|200l [20073) 
searched the Sloan Digital Sky Survey (SDSS) archives to 
find AGNs with low-mass black holes and have found nu- 
merous examples of Seyfert 1 galaxies with A/bh < 10^ 
Mq, which we will refer to as the GH sample. How- 
ever, most of the SDSS objects found by Greene & Ho 
are much more distant than NGC 4395, making detailed 
studies of their nuclei and host galaxies more difficult. 
Although the GH sample does contai n examples of AGNs 
in low-mass disk galaxies (see also iDong et all l2007t ) , 
NGC 4395 remains an exceptional case due to its ex- 
treme late-type morphology as well as its proximity. 

POX 52, also known as G1200-2038 or PGC 038055 
{D = 93 Mpc for Hq = 70 km s^^ Mpc"^), is an inter- 
esting counterpart to NGC 4395. It was originally discov- 
ered in an objective prism study bv lKunth et al.l (|1987f l. 
who classified it as a dwarf disk galaxy with a Seyfert 
2 nucleus, but noted the presence of a weak broad com- 
ponent of the H/3 line. Barth et al. (2004) re-observed 
POX 52, classifying it as a type 1 AGN based on clear 
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Fig. 1.— HST ACS/HRC images of POX 52 and GALFIT models. The top and bottom rows show the F814W and F435W images, 
respectively. In each row, the left panel is the ACS image, the middle panel is the full GALFIT model, and the right panel shows the fit 
residuals. For the F435W filter, the double-Sersic host galaxy model is shown. 



detections of broad Ha and H/3 and noting that its op- 
tical emission-line spectrum is nearly identical to that of 
NGC 4395. Both NGC 4395 and POX 52 can be classified 
as narrow-line Seyfert 1 (NLSl) galaxies, since they sat- 
isfy the defining criterion of having FWHM(H/3)< 2000 
km s~^. NLSl galaxies are often, but not always, char- 
acterized by strong Fe II emission and relatively weak 
[O III] emission when compared to H/3 along with steep 
soft (0.5 — 2.0 keV) spectral slopes and rapid variabihty i n 
the X-ray (jOsterbrock fc Poggdll985HBoller et al]ll996f ). 
NGC 4395 and POX 52 are therefore somewhat atypi- 
cal NLSl galaxies in that they have very high-excitation 
spectra with large [O III]/H/? ratios and very weak Fe II 
emission. 

Using th e vir ial scaling relation calibrated by 
iKaspi et al.l ()2000[ ). Barth et al. estimated a black hole 
mass of Mbh ~ 1.6 x 10^ M©, making this one of the least 
massive black holes identified in any AGN. This mass is 
consistent with the mass expected based on the Mbh — 
relation (extrapolated to small cr») for the measured stel- 
lar velocity dispersion cr, = 36 ± 5 km s~^. They also 
estimated the bolometric luminosity (Lboi ~ 2 x lO"*"^ 
ergs s""'^) and the Eddington ratio (iboi/-^Edd ~ 0.5 — 1) 
for POX 52. Thus, although POX 52 and NGC 4395 
have very similar black hole masses and emission-line 
spectra that are nearly identical in appearance, the AGN 
in POX 52 is ~ 300 times more luminous than that in 
NGC 4395, which ha s an estimate d iboi/^Edd ~ 10^^ 
(jPeterson et al.ll2005D . iBarth et all (2004 ) used ground- 
based images to study the morphology of POX 52 and 
found that the host galaxy was best fit with a Sersic 
model with an index oi n — 3.6 ± 0.2, effective radius of 
— 1'.'2, and absolute magnitude My = —17.6. Com- 



bined with (T, measured from the spectrum, these results 
show POX 52 to fall near the dwarf elliptical sequence 
(Geh a et al ll2003l ) in th e fundamental plane projections 
of iBurstein et al.) (|1997| ). making it the first example of 
a dwarf elliptical to host a Seyfert nucleus. Although 
the past few years have seen a large increase in the num- 
ber of known AGNs with low-mass black holes, POX 52 
remains a valuable target for further study, because of 
its relatively small distance and the fact that it is a less 
luminous and presumably less massive galaxy than the 
majority of the host galaxies in the GH SDSS sample. 

Among nearby dwarf ellipticals, there is little evidence 
for the presence of central black holes. The closest 
well-studied example is NGC 205 in the Local Group, 
and stellar-dynamical observations with the Hubble Space 
Telescope (HST) have set an u pper limit to the bla ck hole 
mass of AfBH < 3.8 x 10"* Mq ijValluri et aI1l2005[ ). which 
is below the mass expected from a simple extrapolation 
of the Mbh — o^^, relation. For dwarf elliptical galaxies 
outside the Local Group, stellar-dynamical observations 
lack the spatial resolution needed to detect black holes, 
and AGN surveys are the best available alternative for 
black hole searches. A double nucleus in one Virgo dwarf 
elliptical galaxy has been interpreted as morphological 
evide nce for a stellar d i sk in orbit around a central black 
hole (jPebattista et al.ll2006f) . but aside from NGC 205 
there are essentially no direct observational constraints 
on the black hole occupation fraction in dwarf elliptical 
galaxies. POX 52 represents a unique opportunity to 
study, in depth, the unambiguous presence of a low-mass 
central black hole in a nearby dwarf elliptical galaxy. 

Here, we present new multiwavelength observations of 
POX 52, including the first targeted HST and X-ray 
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observations of this unusual object. HST images are 
used to examine the host galaxy morphology; our re- 
sults confirm the previous conclusions that the galaxy 
is consistent with a dwarf elliptical morphology with a 
Sersic index of ~ 4. We describe X-ray observations ob- 
tained with both the Chandra and XMM-Newton X-ray 
Observatories in order to study the temporal and spec- 
tral properties of the active nucleus, and to constrain its 
bolometric luminosity and black hole mass. We also in- 
clude observations from the Very Large Array (VLA) and 
^pzteer Infrared Spectrometer (IRS), as well as archival 
data from 2MASS and GALEX, in order to investigate 
the spectral energy distribution (SED) of POX 52 over 
a wide range in frequency. We use recently updated ver- 
sions of the broad-line virial scaling relations to obtain 
new estimates of the mass of the black hole, and in combi- 
nation with the multiwavelength SED we determine the 
bolometric luminosity and Eddington ratio for the AGN. 
In addition to the HST data presented here, we were 
also awarded time as part of the same project to obtain 
ultraviolet and optical spectra of POX 52 with the Space 
Telescope Imaging Spectrograph (STIS), but STIS failed 
before the observations were scheduled to be taken. 

2. OPTICAL OBSERVATIONS 

Observations of POX 52 were made on 2004 Novem- 
ber 18 UT with HST using the High Resolution Camera 
(HRC) on ACS. Four images of 590 s each were taken 
through the F814W (/-band) filter and four images of 
622 s each were taken through the F435W (5-band) fil- 
ter. These images were dithered using a 4-point box pat- 
tern with 0'.'15 spacing to enable the removal of cosmic 
ray hits and bad pixels. An additional 100 s image was 
taken in each filter in the event that the nucleus satu- 
rated in the longer exposures, though this saturation did 
not occur in any image. The images for each filter were 
combined using the Multidrizzle task in PyRAF (version 
1.3) to remove cosmic ray hits and to correct for the ACS 
geometric distortion, and the output images were not sky 
subtracted or rotated during this processing. These HST 
images (Figure [Ij show POX 52 as a bright point source 
centered in a smooth host galaxy with no obvious signs 
of spiral structure in cither filter. 

The two-dimensional fitting program GALFIT (version 
2. Id, Peng et al. 2002) was used to model the host galaxy 
morphology. The central point source in the center of 
the host galaxy was modeled by a point spread func- 
tion (PSF) image, and the host galaxy itself was mod- 
eled by a Sersic component, or combination of Sersic 
components. Since the sky background was not sub- 
tracted in the image processing, a sky component was 
also used within each galaxy model. All magnitudes are 
in the Vegamag system calculated using the zeropoints 
calibrated for HRC^ which are 25.188 for F435W and 
24.861 for F814W, and corrected for Galactic extinction 
using the IRAF package synphot for E{B — V) — 0.052 
mag (jSchlegel et al.lll998[) . 

2.1. PSF Modeling 

The PSF image used in the GALFIT modeling can 
either be created synthetically from a program such as 

* http:/ /www. stsci.edu/hst/acs/analysis/zcropoints 
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Fig. 2. — Radial profile plots of the ACS/HRC standard star 
images along the semi- major axis (SMA), with the Tiny Tim syn- 
thetic PSFs over-plotted in red. The magnitude scale is normalized 
to zero mag arcsec"^ in the PSF core (except for the 2006 obser- 
vation of GD71 which uses the same offset as the 2004 observation 
in order to demonstrate time variability of the PSF structure). 
The bump at ~ 0"4 in the F814W radial profiles is the the long- 
wavelength scattering halo in ACS/HRC. 

Tiny Tim (lKrist|[l993l ) or cropped from a data image us- 
ing the same detector and filter set as the object being 
modeled. Both methods were tested in order to deter- 
mine the most appropriate PSF for our analysis. We 
began by creating a synthetic PSF for each filter us- 
ing Tiny Tim. We centered the PSF in the same po- 
sition as the point source in our galaxy images in order 
to properly reproduce the geometric changes involved in 
the Multidrizzle process. These PSF images were then 
Multidrizzled in the same manner as the object data, 
since the Tiny Tim PSFs created for HRC include the 
geometric distortions in raw HRC images, and therefore 
must be subject to the same corrections to undistort the 
images. In order to evaluate how well the Tiny Tim 
PSF compares to the PSF of a star observed with HRC, 
we searched the archive for unsaturated images of HST 
standard stars observed in the F814W and F435W fil- 
ters, finding three that matched our criteria: G191-B2B, 
GD153, and GD71. Table [1] lists the observation dates 
and exposure times of these three white dwarf stars, in 
which they were located near the center of the CCD. 
The specific observations were chosen to bracket the POX 
52 observation in order to minimize any changes to the 
PSF over time, with the average observation occurring 6 
months from that of POX 52. All of the PSF star images 
were processed through Multidrizzle in the same manner 
as our POX 52 data with the exception that the standard 
star images were sky-subtracted during the Multidrizzle 
process in order to be used as input PSF images with 
GALFIT. 

In order to determine how well Tiny Tim models the 
stellar PSF, we used the Tiny Tim model as the input 
PSF in GALFIT and fit each of the standard stars in both 
of the observed filters. Not only can we examine how 
the shape of the PSF changes for different stars (due to 
color and/or CCD location), but the existence of multi- 
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TABLE 1 
Stellar PSF Observations 



Star 


V (mag) 


Obs. Date (UT) 


Exposure Time (s) 








F435W 


F814W 


G191-B2B 


11.78 


2003 Aug 30 


3.2 


4.0 






2003 Nov 8 


4.0 


14.0 






2005 Feb 28 


4.0 


10.0 


GD153 


13.40 


2004 Feb 15 


14.4 


40.0 






2005 May 19 


14.4 


30.0 


GD71 


13.06 


2004 Feb 3 


12.0 


35.0 






2006 Nov 20 


11.0 


26.0 



pie observations of each of these stars in the HST archive 
allows us to determine the extent of changes to the PSF 
over time. Figure [2] shows radial profile plots made us- 
ing the IRAF task ellipse for the star images compared 
to the GALFIT model using the Tiny Tim PSFs. The 
three left panels are profiles from the three different stars 
found in the HST archive. One can see how the shape of 
the PSF's central peak can vary depending on the star 
used. The two observations of GD71 were taken about 
a year and a half apart. There is not a large change in 
the PSF shape during this time period, but small differ- 
ences are apparent. The differences in profile shapes are 
more apparent in the F435W filter, showing Tiny Tim 
has more difficulty modeling a PSF in this filter. The 
Tiny Tim F435W PSF tends to have a flatter inner core, 
leaving a ^ 0.3 mag arcsec^^ residual out to ~ 0.05". In 
the ACS/HRC camera about 10 - 20% of the PSF flux 
at long wavelengths (A > 7500 A) scatters into a halo 
surrounding the PSF core making PSF mode ling more 
difficult in red passbands fSiria nni et alU lQQS*). such as 
the F814W filter. The halo radius increases with wave- 
length, and Tiny Tim attempts to model its wavelength- 
dependent structure. Despite the red halo (seen as the 
bump at ~ 0.3" in Figure [5]) and the variations seen be- 
tween the stellar PSFs, Tiny Tim did a reasonable job 
modeling each of the standard stars in the F814W filter, 
with few significant deviations between the shape of the 
Tiny Tim PSF and that of the observed star. Since the 
Tiny Tim PSFs reasonably fit the central regions of the 
observed stars in both filters and have the advantage of 
being noise-free, we selected the corresponding Tiny Tim 
PSF to use in the GALFIT modeUng of POX 52. 

2.2. AGN and Host Galaxy Modeling 

We fit the F814W HRC image of POX 52 over a ~ 169 
arcsec^ (491 x 491 pixel) region with a simple model 
consisting of an AGN point source and a single Sersic 
component, where all parameters (including position and 
magnitude for both the point source and Sersic compo- 
nents and effective radius, Sersic index, axis ratio, and 
position angle for the Sersic component) were allowed 
to vary. Figure [3] shows a radial profile plot of the best 
fit from GALFIT with a Sersic index of 4.3 and a re- 
duced of 15.0 for 241,068 degrees of freedom (dof). 
The top panels of Figure [T] show the F814W HST im- 
age (left panel) and the GALFIT model (middle panel) 
with mAGN = 17.5, rrihost — 16.2 and an effective ra- 
dius of 1"08, or 490 pc. The residuals after subtracting 
the model from the galaxy image are shown in the far 
right panel. The strong residuals in the center of the 
image are mostly confined to the PSF-dominated region 
at r < O'.'l and are likely to be the result of PSF mis- 
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Fig. 3. — Radial profile plot of the surface brightness (top panel) 
of POX 52 in the F814W filter with residuals (bottom panel) be- 
tween the data and the total fit. The black points are the data, 
the solid magenta line is the total fit from GALFIT, the dotted red 
line is the PSF component, and the dashed blue line is the Sersic 
component. 

match, but the weaker residuals seen further out from 
the PSF-dominated region of the galaxy could indicate 
small-scale structure in the host galaxy. We added a 
second Sersic component to the model in an attempt to 
minimize the residuals, but it made little difference in 
the overall fit. The additional Sersic component closely 
followed that of the PSF component, leaving the original 
Sersic component to trace the host galaxy as before. The 
relative gain in the overall fit did not justify the need for 
a two Sersic component model, and, therefore, we select 
the PSF-I-Sersic model as the best fit the F814W images. 
Parameters for this fit and the F435W fits can be seen 
in Tabled 

The F435W data proved to be more difficult to model. 
The galaxy image was initially fit using the simple 
PSF-|-Sersic model as above (Figured! l^ft panel), but 
this resulted in an unexpectedly high Sersic index of 
n = 20 (the default maximum value for the Sersic in- 
dex used by GALFIT) and a poor fit compared with the 
F814W model results. In order to get a more reasonable 
result, we added a second Sersic component to the model, 
allowing the parameters for both components to vary 
freely. In the resulting model fit, the first Sersic com- 
ponent has a high central surface brightness that traces 
the central region of the galaxy out to ^ 0"25, resem- 
bling the PSF component but with a slightly broader 
profile, while the second Sersic component has a much 
lower surface brightness that models the majority of the 
host galaxy from ^ 0"5 outward. This outermost Sersic 
component has an index of n = 4.4, similar to the result 
of the F814W fit. 

The effective radius of the host galaxy, as de- 
scribed by the outer Sersic component in the best-fit 
(PSF-f-Sersic-f Sersic) model, is slightly larger than in the 
F814W filter with = 1.22" = 550 pc. The reduced 
is also larger than in the F814W filter (26.2 compared 
to 15.0, even including the second Sersic component in 
the F435W model), and the PSF component has an un- 
expectedly low central surface brightness, as seen in the 
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Fig. 4. — Radial profile plots of the surface brightness {top pan- 
els) of POX 52 in the F435W filter with residuals {bottom panels) 
between the data and the total fit. The black points are the data, 
the solid magenta lines are the total fits from GALFIT, the dotted 
red lines are the PSF components, and the dashed blue lines are 
Sersic components. Left Panel: Best-fit GALFIT model using only 
one Sersic component. Right Panel: Best-fit GALFIT model using 
two Sersic components. 



TABLE 2 
GALFIT Results 



Band 




Component 


Mag 


n 


re (arcsec/pc) 


F814W 


15.0 


PSF 


17.5 










Sersic 


16.2 


4.3 


1.08/490 


F435W 


31.6 


PSF 


19.0 






- Single Sersic 




Sersic 


16.8 


20.0 


5.87/2700 


F435W 


26.2 


PSF 


19.2 






- Double Sersic 




Sersic 1 


18.4 


1.4 


0.060/27 






Sersic 2 


18.2 


4.4 


1.22/550 



radial profile plot (Figure [H right panel). 

The inner Sersic component consists of some combina- 
tion of AGN and host galaxy light, and it is difficult to 
disentangle the relative contributions. The magnitude 
of the PSF component describes a lower limit to the 
magnitude of the AGN of ttiagn = 19.2. We estimate 
an upper limit by repeatedly fitting the galaxy with a 
PSF+Sersic model with a fixed PSF magnitude, increas- 
ing the magnitude each time until it was bright enough 
to account for the majority of the light in the very center 
of the galaxy, while allowing the single Sersic component 
to vary freely as before. This results in an AGN magni- 
tude of mAGN = 18.6, which puts the AGN magnitude 
in the F435W filter in the range of mAGN = 18.6 - 19.2 
mag. We adopt the magnitude of the outermost Sersic 
component as a lower limit to the host galaxy bright- 
ness, mhost — 18.2 in the F435W filter and combine the 
magnitudes of both Sersic components to determine the 
upper limit, mhost = 17.5. This results in a host galaxy 
magnitude range of mhost — 17.5 — 18.2 mag. 

In order to further confirm our choice of PSF, we refit 
the galaxy images with both the PSF-|-Sersic model and 
the PSF+Sersic+Scrsic model using each of the stellar 
PSFs. In each of the fits, regardless of the star used or the 
filter the image was taken in, the PSF component dras- 




FlG. 5.— The Chandra image of POX 52 with intensity (68%, 
90%, and 99%) contours overlaid. Image is 11" on a side. 



tically under-represented the central point source light 
contribution. This verified that the fits in which the 
noise-free Tiny Tim PSF was used produced the best 
overall results. 

Our GALFIT results show that the host galaxy of 
POX 52 has a Sersic index n « 4 in both the F814W 
and F435W filters making it more centrally concentrated 
than a typical dwarf ell iptical galaxy with n « 0.3 — 3 
(iBinggeli fc Jerieiill 19981 ). In a st udy of host ga l axies con- 
taining similar low-mass AGN, iGreene et al.l (|2008f) se- 
lected 19 galaxies from the GH sample for HST imaging 
in order to better determine their host galaxy morphol- 
ogy. From this HST imaging, they find that 60% of the 
sample are disk-dominated galaxies, while the remain- 
der are compact systems with host morphology similar 
to POX 52. These compact galaxies have steep Sersic in- 
dices in t he range of 2 — 4 and reside near dwarf ellip ticals 
from thelGeha et al . (2002), iGavazzi et all (|2005f ). and 
iMatkovic fc Guzman (2005i ) samples in the fundamental 
plane. Although many of the galaxies in the Greene et 
al. HST sample are quite similar to POX 52 in terms of 
Sersic index and overall host galaxy morphology, POX 52 
(^host = -18.4 in the F435W filter) remains less lumi- 
nous than the average g alaxy in the GH sample by '~ 2 
mag ()Greene et al.ll2008l) . 

3. X-RAY OBSERVATIONS 
3.1. Chandra 

We first describe the Chandra observation since the 
source exhibited less spectral complexity in the Chan- 
dra data than in the earlier XMM-Newton observation. 
X-ray observations of POX 52 were carried out using 
the Advanced CCD Imaging Spectrometer (ACIS) in- 
strument on the Chandra X-ray Observatory for a to- 
tal of 25 ks on 2006 April 18 UT. The level 2 event file 
was restricted to the 0.5 — 10.0 keV energy range, and 
there were no periods of high background found. POX 52 
was observed to be consistent with a point source object, 
such that the FWHM of the PSF was - 0'.'5, the 70% and 
90% encircled energy radii were ~ 1.5 and 2 pixels, re- 
spectively^, and the source showed no extended emission 
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Fig. 6. — Left Panel: Light curves of POX 52 taken from Chandra data. The top panel is the full (0.5-10.0 keV) light curve in counts 
s~^, the second panel is the soft (0.5—1.0 keV) X-ray light curve, the third panel is the hard (1.0-10.0 keV) X-ray light curve, and the 
bottom panel is the hardness ratio (HR). Right Panel: Plot of hardness ratio versus the count rate in the full band. 
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Fig. 7. — The Chandra spectrum from 0.5—10.0 keV modeled {thick line) with an absorbed power law (dotted line) and disk blackbody 
(dashed line). Residuals are calculated as Residual = (Data - Model)/Modcl. 



(Figure [5]). Source events were extracted using a cir- 
cular region of radius 10 pixels (4'.'9) centered on POX 
52. The source was too close to the edge of the CCD 
to use an annulus as the background region, therefore, 
a circular background region was extracted 46" away on 
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For ni 



with Chandra 
information, 



PSF 



|http://cxc.harvard.edu/proposer/POG/html/ACIS.html#tth_sEc6.6| 



the same chip with a radius of 24'.'6. There were approx- 
imately 1.2 X 10~^ counts s~^ pixel"^ within this back- 
ground region, giving a background-subtracted average 
source count rate of 0.20 counts s^^. Due to the low 
source count rate, pile-up was not an issue. Data were 
analyzed using standard tasks, such as dmextract and 
psextract in CIAO (version 3.4/4.0/31; CALDB, version 
3.4.0) and spectral models in XSPEC (version 12.3. Ot). 
All light curves were created using time bins of 500 s 
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TABLE 3 
X-RAY Spectral Properties 



Chandra 



XMM-Newton 



C statistic 
Photon Index F 
Disk Tbb (keV) 
Foreground Nu (10^^ 



183 (210 PHA bins, 205 dof) 
1.8 ±0.15 

Q 10+0.03 
'-'■^'-'-0.02 
50+^1 



Partial Covering Nu (10 cm ^) 
Partial Covering Fraction (%) 
F(0.5-10 keV) (lO'^^ ergs cm'^ s" 
L(0.5-10 keV) (10''^ ergs s'^) 



12.6 
1.50 



+0.05 

2.7 
+0.01 

0.37 



111 (104 PHA bins, 97 dof) 

^ ^+0.26 

0.09! 
1.61^ 
58l 
94^ 
3.2^ 
1.38"* 



-0.24 
3+0.03 
-0.01 
1+2.6 
-1.2 
s+8.4 
-9.2 
1+2.2 
-3.3 
,+0.42 
-1.51 
5+0.19 
-0.82 



Note. — All model parameters were fit over a 0.5-10.0 keV energy range and we present observed flux and 
unabsorbed luminosity for both observations. For XMM-Newton, the spectra from both the EPIC-pn and MOS 
cameras were fit simultaneously. 



and split into 0.5 - 10.0, 0.5 - 1.0, and 1.0 - 10.0 keV 
energy ranges. These are denoted as the "full," "soft," 
and "hard" bands, respectively. The full-band light curve 
(FigurelHl top panel) shows strong (factor of 2) variability 
over short (~ 500 s) timescales. Also shown in Figure IH] 
are the light curves for the soft band, hard band, and the 
hardness ratio, HR = (hard — soft)/(hard -I- soft). There 
seems to be no significant correlation between hardness 
ratio and the flux from the full band as seen in the right 
panel of Figure [5] and is confirmed by the results of a 
Spearman rank correlation test. The correlation coeffi- 
cient was —0.16, which describes a weak anti-correlation, 
with a significance of 0.27 that is essentially insignificant. 

The ACIS spectrum was first fit over the range 0.5 — 
10.0 keV with a simple absorbed power law, allowi ng all 
parameters to vary freely and optimized using the ICashI 
(|1979( l C statistic rather than x^- The C statistic is a 
better goodness-of-fit parameter when fitting low signal- 
to-noise ratio spectra. The absorbed power law poorly 
fit the soft end of the observed spectrum, requiring an- 
other component to help model the excess of soft energies 
properly. We added a disk blackbody component to the 
model (PHABS* [POWERLAW + diskbb]), which resulted 
in an acceptable fit to the observed x-ray spectrum (Fig- 
ure [7|) . The best-fitting model had a photon index of 
r = 1.8 ±0.15, and a disk blackbody temperature of 
Tbb — 0.13lg;Q2 keV. The absorbing column density was 
Nh — 0.501J;q;^ X 10^^ cm~^, only sHghtly in excess of the 
predicted foreground column of Nfj ~ 0.4 x 10^^ cm~^ 
based on Galactic H I maps (jPickev h Lockmanl 119901 : 
iKalberla et"ani2005[ ). as determined by the HEASARC 
online Nh calculator-'^^. 

In place of the thermal disk blackbody, we tested 
fitting the spectrum with other thermal components 
plus the absorbed power law component, specifically 
the the rmal bremsstrahlung mod el and the Raymond- 
Smith (jRavmond fc SmithI Il977f) plasma model. The 
bremsstrahlung model, with all parameters allowed to 
vary freely, produced a poor fit that strongly under- 
predicted the soft X-ray end of the spectrum by ~ 3 or- 
ders of magnitude. The Raymond-Smith model produced 
a reasonable fit to the data, similar to the fit when the 
disk blackbody was used, with a C statistic of 182, nearly 
identical to the result for the combined absorbed power 
law and disk blackbody model (C = 183). However, 

-'^'^ http:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi 



this model uses 4 free parameters (3 when the redshift is 
frozen) to fit the data compared to the 2 parameters used 
in disk blackbody model. Due to this difference, we select 
the combined absorbed power law and disk blackbody as 
our best-fit model for simplicity. 

There is possible evidence for an Fe K emission line 
seen in a single bin of energy E = 6.6 keV in the residuals 
of our best-fit model. We added a Gaussian component 
to our model in order to determine constraints on the line 
energy and width. Allowing both the line energy and line 
width to vary freely, but fixing all other model parame- 
ters to our previous best-fit values, XSPEC modeled the 
line centroid at 6.7llg'52 keV with a line width ranging 
from (Tiino = — 0.82 keV and an equivalent width (EW) 
of 0.31q3 keV. The large EW may be explained by the 
X-ray Baldwin effect (see Page et al. 2004, Jiang et al. 
2007), but we cannot draw any fi rm conclusions , since 
the detection is marginal, at best. iNandra et al.l ()2007f ) 
analyzed XMM-Newton spectra of 26 Seyfert galaxies 
(mostly Seyfert 1) and compiled measurements of the Fe 
Ka line. In this sample, the Fe Ka energies ranged from 
5.28 keV to 6.78 keV, with a range of equivalent widths 
of 0.009 — 0.394 keV. Both the measured line energy and 
EW of the Fe Ka from POX 52 are within the range of 
values measured by Nandra et al., although the EW is at 
the high end of the values sampled. The large errors on 
these measurements make this a weak detection at best, 
and deeper observations are needed to obtain a definitive 
detection of the line. 

3.2. XMM-Newton 

XMM-Newton observed POX 52 on 2005 July 8-9 UT 
for a total exposure time of 100 ks on each of the 3 EPIC 
instruments: pn, MOSl, and M0S2. All data were re- 
duced using the Science Analysis System (SAS, version 
7.0.0) and XSPEC following the guidefines of the SAS 
Cookbook and SAS ABC-Guide. For all EPIC instru- 
ments, the source was extracted in a circular region with 
a radius of 32", and due to the proximity of the source 
to the edge of the CCD in the EPIC-pn camera, the 
background was extracted from a circular region free of 
sources, located 107" away from the source with a radius 
of 44"9. Only events corresponding to patterns — 4 (sin- 
gle or double pixel events) were used for the pn event file 
and patterns 0—12 (single, double, triple and quadru- 
ple pixel events) for the MOS event files. All event files 
were restricted to an energy range of 0.5 — 10.0 keV to 
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Fig. 8.— Left Panel: The XMM-Newton EPIC-pn lightcurve (solid line) from 0.5-10.0 keV over the full 100 ks observation. The 
background spectrum is the dotted line. Note the large background flares at 0-10 ks and 40—55 ks. Right Panel: The background- 
subtracted lightcurve of the second half of the observation starting at 57.5 ks, used in the spectral fit. 
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Fig. 9. — XMM-Newton spectra over 0.5-10.0 keV, modeled with an absorbed power law, disk blackbody and a partial covering model 
with all parameters allowed to vary freely. The pn data is in black, MOS 1 is in blue, and MOS 2 is in red. The corresponding models are 
shown as solid lines in their respective colors. The bottom panel shows the residuals (calculated as Residual = [Data - Model] /Model) for 
the EPIC-pn data. 
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better compare with the Chandra data. For aU three in- 
struments, events were excluded that occurred next to 
the edges of the CCD or next to bad pixels. Again, the 
source count rate was low enough to neglect effects due 
to pile-up. 

Again using 500 s time bins, a hght curve (Figured)) 
was created using the entirety of the pn dataset. How- 
ever, the first half of the observation contained two large 
background flares and due to this contamination, only 
the second half of the observation (Figure [51 right panel) 
was used in the spectral analysis (where the total counts 
in the background region are ~ 11.5% of the source 
counts). Examination of the second half of the obser- 
vation shows a similar amplitude of variability to the 
Chandra observations, but over slightly longer (10'^ — lO** 
s) timescales. 

Observations of POX 52 with XMM-Newton yielded 
a significantly lower flux than the subsequent Chandra 
data. Figure [9] shows the spectra from both the pn and 
MOS cameras. There is a hint of a detection of a narrow 
Fe K emission line near 6.4 keV (contained in the bin 
centered at 6.375 keV with a width of 0.125 keV), as 
can be seen in the EPIC-pn residuals (Figure [H lower 
panel). A simple absorbed power law and a combined 
absorbed power law plus disk blackbody both proved to 
be poor fits, even when the same parameter values were 
used from the Chandra fit, under-predicting the spectral 
shape by 1 — 2 orders of magnitude below ^ 1 keV. 

The XMM-Newton spectrum is harder than in the 
Chandra observation and requires more absorption to fit 
the spectral shape. Adding a partial covering component 
to the model (phabs*pcfabs*[powerlaw -|- diskbb]) 
improved the fit significantly. We fit the data using this 
model, allowing all parameters to vary freely. The pho- 
ton index and blackbody temperature of this model are 
consistent within the errors with values derived from the 
Chandra data. The major difference is the need for much 
higher levels of absorption. The absorbing column den- 
sity is three times that needed to model the Chandra data 
and the model adds a partial covering absorbing column 
density of 58t^;2 x lO^^ cm with a partial covering 
fraction of 94^3 3%. Since the partial covering fraction is 
very high, we tested forcing the fraction to 100% or sim- 
ply using one larger absorbing column density, but both 
scenarios tended to over-absorb the power law in the soft 
X-ray end. The unabsorbed luminosity calculated for the 
XMM-Newton observation is within the la uncertainty of 
the Chandra unabsorbed luminosity, which suggests that 
the properties of the central engine itself did not change 
significantly between these two observations. The factor 
of 4 flux difference between the two X-ray observations 
is then explained by the additional absorption needed to 
fit the XMM-Newton spectra. 

4. ULTRAVIOLET OBSERVATIONS 

4.1. XMM-Newton Optical Monitor 

The Optical Monitor (OM) on XMM-Newton took a 
series of 72 images simultaneously with the X-ray ob- 
servation, using the UVM2 filter, which is centered at 
2310 A, extending from ~ 1900 - 2950 A (FWHM of the 
PSF = I'.'S ~ 4 pixels). Each of these images had an 
exposure time of 1 ks. The images were processed using 
the XMM-Newton SAS routine omichain to output cal- 



ibrated photometric data via aperture photometry. Out 
of 72 images, POX 52 was easily detected in 60 images 
by the source detection algorithm in omichain and can 
be seen as the central object in Figure [TOl Five addi- 
tional images were reprocessed with omprep, omdetect, 
ommag, ommodmap, and ommat separately in order de- 
tect POX 52 by varying the detection threshold. For the 
remaining 7 images, the omichain routine was unable to 
detect or measure the source count rate due to high back- 
ground count rates, and these images were discarded. 

We used the photometric results of omichain to con- 
struct a UV light curve, shown in Figure [TlJ The mea- 
sured flux was corrected for Galactic extinction by using 
the [SchlcRcl ct al. (1998) reddening maps and averaging 
the ICardelli et al.l f|1989l) extinction law over the OM 
passband resulting in Aqm = 0.32 mag for the UVM2 
filter. The mean flux density over the set of observations 
is /^(2310 A) = (1.5 ±0.1) X 10-27 ergs cm'^ s'^ Hz"!. 
Inspection of the light curve shows that the source is not 
highly variable in the UV. Only 25% of the individual 
points differ from the mean flux level by more than la, 
and only 3/64 (or 5%) of the points differ from the mean 
flux by more than 2a; there are no Str outliers. The ob- 
served variability is therefore generally consistent with 
the expected level of variations from photon-counting 
statistics for a constant source. To further quantify any 
possible variability, we calculated the normalized excess 
variance o f the light curve, follow ing the prescription de- 
scribed bv lNandra et all ()1997af ): 

1 ^ 

where Xi denotes the count rate of the i-th point in the 
light curve, ai is its uncertainty, /i is the mean of the 
Xi values over the entire light curve, and N is the num- 
ber of points in the light curve. The resulting excess 
variance is very small and actually has a negative value: 
'''nxs = —0.0029. This result confirms that there is no 
evidence for significant variability in the UV light curve. 
The negative value of a^^^ implies the error bars on the 
measurements are most likely slightly overestimated by 
the SAS processing and detection algorithm. 

4.2. GALEX 

As part of the CALEX All-sky Imaging Survey (AIS), 
POX 52 was imaged simultaneously in the near-UV 
(NUV) and far-UV (FUV) bands on 2006 February 9 
UT for a total of 114 s. The NUV and FUV bands span 
the ranges 1771 - 2831 A and 1344 - 1786 A, respec- 
tively, with effective wavelengths of 2271 A and 1528 
A. The FWHM of the PSF for the NUV filter is 4'.'0 
(-2-3 pixels) and for the FUV filter is 5'.'6 (-3-4 
pixels). Figure [10] shows both NUV and FUV images 
of POX 52 alongside one taken with the OM on XMM- 
Newton. Galactic extinction was calculated using the 
same method as for the OM, resulting in Anuv = 0.43 
mag and Apuv = 0.84 mag. These were then used to 
correct the fluxes of POX 52 taken from the AIS cat- 
alog, resulting in a NUV flux density of /j.(2271 A) = 
(1.88 ± 0.08) X 10-27 ergs cm^^ Hz^^ and a FUV 
flux density of ^(1528 A) = (1.82 ± 0.10) x 10^27 
ergs cm-2 Hz-i. The NUV band of GALEX sur- 
rounds the UVM2 filter on the OM, giving us the oppor- 
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Fig. 10. — Ultraviolet images of POX 52 (central object in each panel), 
image taken from our sample of 72 images. 
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Fig. 11. — UV light curve of POX 52 from the Optical Monitor 
onboard XMM-Newton. The horizontal line represents the mean 
observed flux. 

tunity to check for long-term variability in the UV. The 
flux densities calculated from both instruments are con- 
sistent within the 2a uncertainty range, meaning we see 
no significant change in UV luminosity in the 7 months 
between the separate observations. 

In order to estimate the UV host galaxy contribution, 
we take the observed range of NUV— F and FUV— col- 
ors of Virgo dwarf elliptical galaxies from iBoselli et al.l 
(|2005l ) and compute possible NUV and FUV flux densi- 
ties based on the pre vious measurement of Vhost = 17.4 
mag for POX 52 from lBarth et all (|2004D . This results in 
a estimated host galaxy magnitude in the range of NUV 
= 20.4-22.4 mag (AB) and FUV = 21.9-24.9 mag (AB), 
which we convert to a host galaxy NUV flux density of 
/^(2271 A) = (2.5-4.0) X lO'^s ergs cm'^ s^^ Hz^i and 
a FUV flux density of /^.(1528 A) = (4.0 - 63) x lO^^o 
ergs cm~^ llz~^ for POX 52. We take the aver- 
age flux density value in each of these bands to cor- 
rect for host galaxy contamination and find corrected 
AGN flux densities of /,.(2271 A) = 1.85 x 10-^7 
ergs cm-2 s'^ Hz"! and ^(1528 A) = 1.80 x lO^^^ 
ergs cm~^ s~^ Hz~^ for the NUV and FUV, respectively. 
Correcting the OM measurement in a similar manner re- 
sults in a negligible change to the flux density. We note 



that, even for the highest estimate of host galaxy UV 
flux, the source is still very AGN dominated in the UV. 

5. RADIO OBSERVATIONS 

POX 52 was observed for one hour on 2004 October 10 
UT w ith the VLA in A configuration (T hompson et al.) 
fT980l) . as part of a larger program to observe low-mass 
active galaxies from Greene & Ho (2004; see Greene et 
al. 2006). Observations were carried out at 4.860 GHz 
(6 cm; C band) with a bandwidth of 50 MHz for each of 
two intermediate frequencies separated by 50 MHz. For 
purposes of phase calibration, observations of the target 
were alternated every 3-5 minutes with observations of 
a bright radio source (PKS J1159-2148, located - 1° 
away). Our overall fiux scale is tied to 3C 48, with an 
assumed flux density of 5.4 Jy, and uncertainties in the 
absolute fluxes of the phase calibrators, estimated to be 
~ 5%, dominate the flux scale uncertainties. 

After removal of corrupted records, flux and phase cali- 
bration were performed within AIPS (iGre isen '20031. We 
then Fourier transformed the observed visibilities to cre- 
ate an image of intensity on the sky, using a pixel scale 
of 0'.'06 and a resulting synthesized beam of AO 0'.'6. 
We find no confusing sources within the primary beam 
(^ 9'), and so, because deconvolution algorithms such 
as CLEAN may change the noise statistics in a source- 
free image, we have calculated the upper limit of 0.078 
mjy directly from the direct Fourier transform. This 
upper limit represents 3x the measured noise limit in 
the uncleaned map of = 0.026 mJy and corresponds to 
a radio power of P < 8.1 x 10^^ W Hz^^ and a radio- 
loudness parameter {R = fscm/ f^^^^j) for POX 52 of 
R < 0.28, where a radio-loud object is usually defined to 
have i? > 10 ([Kellermann et al. 'l989). 

These values are consistent with the majority of th e 
other GH sample sources studied bv lGreene et al.l (|2006f) . 
Out of 19 objects studied, Greene et al. only signif- 
icantly detected one object, with 3(t upper limits to 
the radio power of the other 18 objects around P < 
a few X 10^^ W Hz^^. The high upper hmits for the 
GH objects relative to POX 52 are due to the larger dis- 
tances to the GH sample. The upper limits on the radio- 
loudness parameter may be a better comparison between 
POX 52 and the GH objects, which range from R < 0.68 
to i? < 9.9 for the GH sample, but the majority of these 
objects are well below the i? > 10 cutoff that usually 
describes a radio-loud object. Stacking the GH sample 
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results in i? < 0.27, giving a stricter upper limit to the 
radio-loudness of the sample as a whole. 

6. INFRARED OBSERVATIONS 
6.1. Spitzer Space Telescope 

POX 52 was observed on 2006 Ju ly 24 UT with the 
Spitzer Infrared Spectrograph (IRS) (jHouck et al.ll2004D 
as part of a sample of low-mass Seyfert galaxies (Pro- 
gram ID: 30119). The analysis of the entire sample, in- 
cluding measurements of emission lines and continuum 
shape, will be presented in a later paper, but we include 
the POX 52 spectrum here in order to supplement our 
multiwavelength coverage of this object. It was observed 
in the Short-Low and Long-Low slit modes, which corre- 
spond to a total wavelength range of 5.2 — 38 fim. Expo- 
sure times were 2 x 60 s each in the SLl and SL2 settings, 
and 2 x 120 s each in the LLl and LL2 settings. Data 
were processed prior to download with IRS Pipeline Pro- 
cessing version 15.3.0. Data were first cleaned using the 
IDL routine IRSCLEAN (v. 1.9), and the images were 
then coadded following the procedure in the Infrared 
Spectrograph Data Handbook^^ (v. 3.1). The spectrum 
from each slit was extracted and calibrated using SPICE 
(v. 2.1.2). For each slit, the two nod positions were sub- 
tracted from each other in order to sky subtract the re- 
sulting data. Spectra from each of the four slits were 
combined and any overlapping data from adjacent wave- 
bands were averaged together in the final product. Data 
from the LL slit was multiplied by a factor of 1.4 in order 
to properly match data from the SL slit. Figure [12] shows 
the combined spectrum with the most prominent emis- 
sion lines identified, including [Ne II] A12.81 /im, [Ne III] 
A15.56 /im, [Ne V] A14.32 ^m, [O IV] A25.89 ^lm, [S III] 
AA18.71, 33.48 ^m, [S IV] A10.51 fim, [Si II] A34.82 ^m, 
and PAH features at 6.2 ^m, 7.7 ^m, and 11.3 fim. 

The IRS spectra of POX 52 shows possible weak sili- 
cate absorption near 10 /im seen in some Seyfert 1 galax- 
ies, but more commo nly associated with Seyfert 2 galax- 
ies (|Hao et al.l l2007t ) . Assuming a power law fit over 
the 30 — 20 /xm range with oc i^^", POX 52 shows 
a steeper spectral slope ctsp 20 ~ 2.1 than many Seyfert 
1 galaxies (Buchana n et al.ll2006D . The spectral slope 
from 15 — 6 /zm of ais-e ~ 1-16 is significantly flatter 
than a3o_2o suggesting the best continuum fit to be a 
bro ken power law, which i s often seen in Seyfert 1 spec- 
tra (jWeedman et al.l[2005h . The presence of PAH fea- 
tures and the observed red continuum (seen as the steep 
Q^30-'20 slope) suggests possible star formation in the host 
galaxy, however, we defer any further analysis to the later 
paper on the full Spitzer sample. 

6.2. 2MASS Point Source Catalog 

Near-infrared magnitt ides were taken from t he 2MASS 
Point Source Catalog (jSkrutskie et al.l l2006f) in the J 
(15.7 mag), H (15.0 mag), and Ks (14.5 mag) bands to 
be included along with our other data in the SED (Fig- 
ure I15p . These magnitudes include flux from the host 
galaxy as well as the AGN. In order to correct for this 
we used averaged B — J = 2.05, B — H = 2.57, and 
B — if = 2.67 color indices derived from photometry o f 
10 dwarf ellipticals in the Virgo Cluster (| James! [19941 ). 

http:/ /ssc. spitzer. caltech.edu/irs/dh/ 




Fig. 12.— Spitzer IRS spectrum of POX 52 from 5-38 ^tm. 
The most prominent spectral lines are identified. The inset shows 
a close-up of the 5 — 12 /xm region. 

Our GALFIT estimate of the i3-band mhost (corrected 
from the F435W ACS/HRC flfier to the standard John- 
son B filter using the IRAF routine calcphot) was used in 
order to estimate the JHK contribution from the host 
galaxy. The host galaxy contributions were then sub- 
tracted from the 2MASS magnitudes, resulting in the 
following AGN point source JHK magnitude estimates: 
Jagn = 16.8l[J:^, -^AGN = 15.81°^, and i^ACN = 
IS.Olg'g mag. The errors were calculated based on the 
range oi B ~ J, B — H, and B — K colo rs for dwarf 
elliptical galaxies derived from I James! (|1994( l. 

7. DISCUSSION 
7.1. Host Galaxy Structure 

iBarth et all (|2004f) previously classified POX 52 as a 
dwarf elliptical based on its morphology, their measure- 
ments of Mb = —16.8, Te = 1'.'2, and n = 3.9, and 
the location of POX 52 in the "K-spac e" projections 
("Ben der et all 11992! : !Burstein et"an!l997D of the funda- 
mental plane. These k parameters combine stellar ve- 
locity dispersion, surface brightness, and effective radius 
together in such a way, which correspond to the physical 
parameters of stellar mass and mass-to-light ratio, that 
a clear split can be seen among giant ellipticals, dwarf 
ellipticals, dwarf spheroidals, and globular clusters. In 
these projections, POX 52 falls nearest th e Virgo dwarf 
elliptical galaxies from iGeha et al.! ([2003), and is offset 
only slightly toward a higher surface brightness and a 
lower mass-to-light ratio. Our new HST observations do 
little to change the location of POX 52 in these projec- 
tions, as the surface brightness and effective radius found 
in this paper are consistent with the previous values. 

The addition of similar objects from the GH sam - 
ple, based on new HST imaging ([Greene "eTalllMiah . 
to these fundamental plane relations suggests that they 
are more centrally concentrated than typical dwarf ellip- 
ticals. However, in plots of vs. My and Tg vs. My, 
they reside in between the distributions of dwarf ellipti- 
cals and classical ellipticals, as a possibly intermediate 
population. These galaxies clearly have properties that 
are in common with both populations, but they follow 
the distribution of dwarf ellipticals when kinematic data 
are combined with photometric data in the "K-space" 
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projections and therefore, we conclude that POX 52 and 
similar objects are most closely related to dwarf ellipti- 
cals, although they remain unusual for having unexpect- 
edly highly Sersic indices. 

We note a difference in terminology b etween what some 
authors call a dwarf ellip tical galaxy (jBinggeh fc Jerieiil 
119981 iGeha et all |2002| ) and a different designation, 
spheroidals, in order to distinguish the typically less cen- 
trally concentrated ob jects from low-lum inosity versions 
of classical ellipticals ()Kormen d^IHsg). In this paper, 
we refer to these galaxies as dwarf ellipticals for consis- 
tency with the nom enclature adopted in our earlier paper 
(jBarth et al.ll2004r ). 

From our GALFIT decompositions, we can conclude 
that the host galaxy is well fit in both of the observed 
filters by a Sersic profile with n « 4 and « 1" at r > 
0'.'5. The structure of the inner region of the host galaxy 
remains complicated due to the second Sersic component 
needed to provide a reasonable fit in the F435W filter. 
We noted earlier that the inner Sersic component seems 
to be compensating mostly for the PSF mismatch, but 
its broader profile suggests the possibility of a spatially 
compact, blue component, such as a nuclear star cluster. 
Looking at the B — I color of POX 52 with and without 
the inner Sersic component included, we can bracket the 
range of possible colors of the host galaxy. Excluding 
the inner Sersic component, we find B — I ~ 2.0 for the 
host galaxy, which is similar to fiducial values (B — I = 
2.0 — 2.1) for so/Sab galaxies found in iFukugita et al.l 
(|1995f ) a nd from synph ot calculations [B — I — 2.2 — 2.3) 
with the iKinnev et al.) (1996) SO/Sab galaxy templates. 
If we assume that the compact component is entirely due 
to starhght, th e resulting color of B — I = 0.74 is bluer 
than both the iFukugita et al.l ()1995l ) fiducial value for 
Im galaxies {B — I = 0.9)and all of the starburst galax y 
templates (B - / = 0.9 - 1.6) from Kinn ev et all (jl996f ). 
Therefore, we conclude that the inner Sersic component 
in the F435W image must contain a substantial amount 
of light from the AGN point source. However, we cannot 
determine the exact fractions of AGN and starlight in 
this component. This leaves a wide range of possible 
colors for the host galaxy. 

7.2. X-ray Spectral Properties 

X-ray spectral variability is not an unusual occurrence 
in an AGN, but is normally manifested in small changes 
in the spectral slope or in the width or amplitude of an 
emission line. Strong variability in the absorbing col- 
umn density is more rare, but is often attributed to par- 
tial covering. In this scenario, there exists a patchy ab- 
sorber surrounding the AGN made up of many smaller 
gas clouds orbiting the central engine with close to Ke- 
plerian velocities. The variations seen in column density 
are therefore related to the number of clouds along the 
line of sight. This has become a common model to use 
when a change in the spectral slope on the soft end is 
observed, as has been see n in several NL SI galaxies, in- 
cludin g IRAS 13224-3809 (IBoller et al.l[200|jGallo etaH 



|2004£]). RX J22 17 9-5941 (IGrupe et al.ll2004aD. Mr k 335 



(|GruT)e et al.ll2007l). Mrk 1239 (iGrupe et al.ll2004bD . and 
IH 0707-495 (IGallo et al.ll200m 

The bright NLSl Mrk 335 has recently experienced a 
change in absorption very similar to what we observe in 
POX 52. Mrk 335 has regularly been studied in the X- 



rays since 1971 with no dramatic changes in flux until 
2007, when IGrupe et ;ill (|2007D measured a factor of 30 
decrease in flux compared to all previous observations. 
When the spectra were modeled, Grupe et al. were un- 
able to obtain a good flt with scaled-down versions of 
previously used models and found that adding a partial 
covering model to an absorbed power law with a high 
partial covering fraction of 0.9 — 1, much like that of 
POX 52, enabled an acceptable flt. When corrected for 
this additional absorption, the change in intrinsic lumi- 
nosity is only a factor of 4 — 6 when compared to earlier 
observations. Grupe et al. note that a partial covering 
model can also reasonably fit their 2006 data with a lower 
partial covering fraction of ~ 0.5, but does not result in 
the best fit for their 2000 data. This suggests the partial 
covering absorber was moving into the line of sight as 
early as 2006. 

Assuming the sizes of these discrete absorbing clouds 
do not vary much from one source to another, and that 
the accretion disks and surrounding X-ray-emitting re- 
gions of less massive black holes like POX 52 are, in fact, 
scaled-down versions of those associated with their more 
massive counterparts, one would expect the X-ray ab- 
sorbing region to be less voluminous and therefore con- 
tain fewer absorbing clouds. Furthermore, if the differ- 
ences in column density are attributed to the number of 
clouds instead of the size of clouds, it is possible that an 
X-ray-emitting region surrounding a black hole the size of 
POX 52 would have a higher likelihood of being observed 
with almost no absorbing clouds in the line of sight than 
for a more massive source due to the differences in the 
number of clouds associated with each source. The black 
hole in POX 52 is ~ 2 ord ers of magnitude les s mas- 
sive than that in Mrk 335 (iPeterson et al.l [19981 12004 
IVestergaard fc Peterson! l2006t ) and moved from an un- 
absorbed state to the partial covering state in less than 
9 months, whereas the change in a bsorbing column i n 
Mrk 335 occurred over ~ 1 — 7 years (IGrupe et al.ll200l . 
These observations give upper limits to the timescales 
of absorption variability based on the partial covering 
s cenario. 

iBarth et al] (|2004D examined the ROSAT all-sky sur- 
vey image of the position of POX 52 with an exposure 
time of ^ 90 s. Although there was no prominent source 
at the position of POX 52, Barth et al. showed there were 
two counts above the background rate within a circular 
aperture of radius equal to ~ FWHM of the PSF of the 
detector, for a 3.0 — 3.5fT detection. With this in mind, 
we used the WebPIMMS^^ online calculator to determine 
what state POX 52 might have been in when this ob- 
servation was taken. Inputting our best-fit models from 
XMM-Newton and Chandra, we find that the unabsorbed 
state would have resulted in a 0.4 - 2.4 keV i?OS'^T de- 
tection in the range of 4 — 14 counts, and the absorbed 
state would have produced 1 — 4 counts during the ob- 
servation. From this we conclude that the ROSAT flux 
is consistent with the XMM-Newton absorbed state, or 
possibly a lower luminosity state. Further observations 
would help determine the more common X-ray state of 
POX 52 and better constrain the timescales of this vari- 
ability. 

We see no change in the intrinsic power law spectral 

^'^ http:/ /heasarc. nasa.gov/Tools/w3pimms. html 
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slope of r « 1.9 between these two observations. This is 
somewhat unusual given the large amount of both short- 
term 500 s) and long-term (~ 10'' s) variability that 
characterizes many NLSl galaxies, including POX 52. 
For example, POX 52's most similar companion, NGC 
4395 has been shown to have strong variations not only 
in flux (factor of 2 over 5 months, O'Neill et al. 2006), 
but also in spectral slope (P w 0.6 — 1.7 over 2 years, 
Moran et al. 2005). A 2 - 10 keV spectral slope of 
0.6 is unusually low for either a Seyfert 1 galaxy or a 
NLSl galaxy; a slope of ~ 2 is far mo re common for 
Seyfert 1, NL SI gal a xies, and quasars (ICrummv et al.l 
120061 iLeighlvl [l999bl : iPiconcelh et al.ll2005l ). The tem- 
perature of the blackbody component, used to model the 
soft excess seen in many NLSl galaxies, also remained 
unchanged for POX 52 between the two observations, 
with a temperature of kT — 0.13 keV that is within 
the normal ran ge (0.1 — 0.3 keV) seen in NLSl galaxies 
()Leighlvlll999bf ). When NLSl spectra with soft excesses 
are fit with a model consisting of only a power law or ab- 
sorbed power law com ponent, the r esult is a steep spec- 
tral slope of r « 3 — 5 (jBoller et al.l[T9 96:), as opposed to 
the more common NLSl spectral slope of P ~ 2. 

In comparison, the similar low-mass GH sample ob- 
jects observed wi th Cha ndra also have slightly different 
X-ray properties ([Green e & Ho 2007a). These 10 objects 
were observed for 5 ks each, resulting in only 5 objects 
with enough counts to allow for spectral fitting, and even 
this could only be done in the 0.3 — 5 keV range. A soft 
photon index was derived from the hardness ratios for 
each object, giving a range of 1 < P5 < 2.7 for the sam- 
ple which brackets the 0.5 — 10.0 keV photon index for 
POX 52. The GH sample objects show no signs of either 
a soft excess or intrinsic absorption. The observations 
did not allow for study of short-term variability due to 
lo w count rates. 

iMiniutti et"al] ([2008') re-observed four of these GH ob- 
jects with XMM-Newton for 20 - 30 ks, in order to obtain 
a high enough signal-to-noise ratio for spectral fitting. 
These objects have remarkably similar spectra in the 
0.5 — 10.0 keV range to the Chandra spectrum of POX 52, 
and fits with power law plus disk blackbody models pro- 
duce photon indices in the range of 1.8 < P < 2.5 and 
temperatures ~ 0.11 keV. Although the blackbody tem- 
perature of the soft excess is similar to what is seen 
in both the Chandra and XMM-Newton observations of 
POX 52, the hard photon index of these GH objects is, 
on average, slightly harder than what has been observed 
in POX 52. These objects are an order of magnitude 
more luminous in the X-ray than POX 52, and contain 
black holes with ma sses around an order or magnitude 
larger than POX 52 (j Greene fc Holl2007b[ ). 

7.2.1. Power Spectral Density 

AGNs are highly variable in X-rays and with good data 
sampling, one can study the frequency dependence of 
these variations by computing the power spectrum den- 
sity (PSD), which is a measure of the amount of variabil- 
ity at a given frequency. The PSDs of AGNs and Seyferts 
are often fit with a combination of two power laws with a 
slope of —2 at the high-frequency end and then breaking 
to a slope of —1 at the low- frequency end. The frequency 
at which this change occurs is known as the break fre- 
quency, i^B , and appears to correlate inversely with black 
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Fig. 13. — PSD derived from Chandra data using a light curve 
with 100 s time bins. The PSD itself has been binned by a factor 
of 5. A power law with slope —0.95 is shown as a solid line. 

hole mass (|Markowitz et al.l[2003[ ). 

Since the Chandra data have a higher signal-to-noise 
ratio than the XMM-Newton data and more closely rep- 
resent the intrinsic state of the central engine without 
the heavy obscuration seen with XMM-Newton, we chose 
to carry out the PSD analysis on the Chandra dataset. 
We start by creating a light curve using 100 s time bins 
covering the 0.5 — 10.0 keV energy range. The Xronos 
(v5.21) tool powspec (vl.O) was used to create the PSD 
from the light curve and can be seen in Figure fT3l binned 
by a factor of 5 in frequency in order to increase the 
signal-to-noisc ratio. There may be a possible flattening 
of the PSD from low frequencies and a break to a steeper 
slope at ~ 3 X 10"'' Hz. However, the PSD is easily fit 
with a single power law with a best-fit slope of —0.94 
and a = 2.87 for 8 dof. This fit is good enough that 
a broken power law fit is not justified. To test whether 
any constraints on vb could be derived, we attempted 
to fit a broken power law to the data. Unfortunately, all 
attempts at fitting the data while allowing parameters to 
vary freely produced either a break frequency outside the 
range of our data or a broken power law with no break, 
and the result is that both power law models had identi- 
cal slopes. Arbitrarily forcing a specific break frequency 
produced fits to the data with low-frequency slopes of 
~ —0.3 and high-frequency slopes of —1.4, neither of 
which approximate the expected slopes of —1 and —2 
seen in most AGNs and Seyferts. From this we conclude 
that if is actually within the range detectable from 
our data, that our PSD is of such low quality that we 
are unable to measure it, or alternatively, vb lies at a 
frequency higher than ~ 5 x 10~^ Hz. 

Even without an accurate measurement of i^b, we can 
still gain information from the relationship between ex- 
cess variance and Lx- There seems to be a noticeable 
split between the distribution of cr^xs foi" NLSl and clas- 
sical Seyfert 1 galaxies, with the excess variance of the 
NLSl galaxies consistently larger th an that of the S eyfert 
1 galaxies at a given luminosity (Lcighlyl Il999al ) . To 
test this for the low-mass case of POX 52, we calcu- 
late the excess variance from the entirety of the usable 
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Fig. 1 4 . — Plo t of o-„yB vs. 2-10 keV luminosity. NLSl data from 



ILeighiy fl999a'!B) arc shown as filled circles. Classical Seyfert 1 
data from Nandra ct al. (1997a b) are shown as open circles. POX 
52 is shown as a red (Chandra data) or a blue (XMM-Newton 
data) cross. The unabsorbed luminosity is used for the XMM- 
Newton data point in order to represent the intrinsic luminosity 
of the source. Similar l ow-mass o bjects from the GH sample are 
shown as open triang les l IMiniutti et al..2008i) . including NGC 4395 
lllwasawa et al.ll2000l) as a filled triangle. 

light curve for both our XMM-Newton (second half of 
the observatfon, ~ 50 ks) and Chandra (~ 25 ks) data, 
using 500 s and 100 s bins, respectively and plot the re- 
sults against the 2 — 10 keV unabsorbed luminosity, along 
with NLSl data fro m Leij^hlvl (|1999allbl ). Seyfert 1 data 
fromlNandra et ahl ( 1997a't?). and GH object data from 
iMiniutti et all ()2008i ) in Figure [Ml From this analysis, 
we find the excess variance of the two observations to 
be quite similar, a'^^^{Chandra) = 0.0388 ± 0.0003 and 
al^^{XMM-Newton) = 0.0494 ± 0.0073, despite the sub- 
stantial additional absorber present in the XMM-Newton 
observation. This suggests that in the case of POX 52, 
the variability seen in the XMM-Newton observation is 
due to intrinsic variability and not due to fluctuations in 
the partial covering or the absorbing column density. 

7.3. Spectral Energy Distribution 

We combine all of our data and present one of the most 
complete SEDs (Figure [151 Table [4]) of a low-mass AGN 
to date, covering the radio to t he X-ray. For com pari- 
son, we include the templates oflElvis et al.l (11994 ) from 
Palomar-Green (PG) quasars f Sch midt k Greed ll983D . 
normalized to the flux of our NUV GALEX data point. 
The data from our observations closely match those of 
the radio-quiet (RQ) template for the majority of the fre- 
quency range shown. The Chandra data show POX 52 as 
X-ray bright compared to the RQ template and we again 
note that the XMM-Newton data represent a period of 
high absorption, not a decrease or change in the flux of 
the central object itself. It should also be noted that the 
observations are non-simultaneous and could therefore 
represent POX 52 in various states of emission, so care 
should be taken in examining the SED and comparing 
fluxes from different observations. 

We approximate the shape of the SED near the OM 
observations as a power law with f^ociy^^ and extract a 
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Fig. 15.— SED of POX 52. The dashed lines are the radio- 
loud (short dash) and radio-quiet (long dash) SED templates from 
lElvis et al.l ((T994) normalized to the UV data point. The X-ray 
data is in green with the dotted (Chandra) and dashed (XMM- 
Newton) lines representing the Icr uncertainty range for the respec- 
tive observatories. The HST (squares) and OM (asterisk) data are 
in blue, the magenta circles are from 2MASS, the purple triangles 
are from GALEX, the red spectrum is from Spitzer IRS (binned by 
a factor of 5 for clarity), and the gold arrow represents the upper 
limit on the VLA data. The top panel shows a close-up of the IR to 
X-ray data and the bottom panel shows the full extent of the data, 
excluding the strongly absorbed XMM-Newton observation. 

UV luminosity at 2500 A of 1.5 x 10^"^ ergs s^^ Hz^^. 
Combining this with the 2 keV unabsorbed luminos- 
ity derived from the simultaneous observations with the 
EPIC cameras on XMM-Newton, we can calculate aox, 
the slope of a hypo thetical power law conn ecting the 
UV to the X-ray 13 ([Tananbaum et al.lll979f) . We find 
aox — —1-3, and this value does not change upon sub- 
stitution of th e Cha ndra 2 keV unabsorbed luminosity. 
IStrateva et al.l (|2005il studied aox in a sample of 228 
optically selected AGNs and found that aox decreased 
with increasin g luminosity. T his correlation has been up- 
dated by Steff en et al.l (|2006( ) to include higher redshift 
objects, resulting in aox = — 0.137Luv + 2.638, where 
LuY is the luminosity of the object measured at 2500 A. 
We plot the data from Strateva et al. along with that of 
POX 52 and s imilar i ntermediate-mass black holes from 
iGreene fc Hoi (|2007aD in Figure [H showing that POX 
52 falls among its intermediate-mass counterparts and 
within the la errors of the extrapolated correlation be- 
tween aox and UV luminosity calibrated for more lumi- 



POX 52 



15 




26 28 30 

log L„, (2500 A) 

Fig. 16. — Plot of QQx vs. UV luminosity. Black dots are lu- 
minous type 1 AGNs from lStrateva et al.l II2005I ). arrows represent 
uppe r limits, the open tria ngles are intermediate-mass black holes 
from lGreene fc Hcf l l2007al '). and the red cross is POX 52 derived 
from simultaneous UV and X-ray data taken with XMM-Newton. 
The solid line is the linear relation from .Steffen et al., (,200& ) and 
the dashed lines represent the la error on that relation. 



TABLE 4 
SED Data 



Rest Wavelength 


log [uL^/ergs s ^] 


Source 


or Energy 




6 cm 


< 36.127 


VLA 


36 fj,m 


42.58 


Spitzer IRS 


30 fim 


42.39 


Spitzer IRS 


20 Atm 


42.22 


Spitzer IRS 


10 fim 


41.75 


Spitzer IRS 


6 fim 


41.70 


Spitzer IRS 


2.12 /xm 


41.85 


2MASS Ks 


1.61 fim 


41.97 


2MASS H 


1.22 ^m 


41.89 


2MASS J 


7966 A 


41.97 


HST F814W 


4257 A 


41.88 


HST F435W 


2260 A 


42.29 


XMM-Newton OM 


2222 A 


42.40 


GALEX NUV 


1495 A 


42.56 


GALEX FUV 


1.0 keV 


41.62 


Chandra 


2.0 keV 


41.59 


Chandra 


10.0 keV 


41.71 


Chandra 



nous type 1 AGNs from StefF en et al. 

Using the lElvis et al.l (| 19941 ) UV bolometric correction, 
we estimate Lboi ~ 1.1 x 10^^ ergs from the OM ob- 
servations. Repeating this with the Elvis et al. S-band 
correction and the HST measurement, results in Xboi 
« 6.5 X 10**^ ergs s^^, which is a factor of ~ 2 smaller 
than that of the UV data. Both of these estimates are 
based on observations of POX 52 in a single band, con- 
stituting a single point each on the SEP. Using the tem- 
plates derived bv lElvis et all (|1994[ ). shown on the SED, 
we estimate the Lboi using all of the data collected. With 
the templates normalized to match the GALEX NUV lu- 
minosity of POX 52, the radio-quiet template matches 

^'^ If a power law of the form fv oc !y°'Ox js assumed, aox is 
then defined by oqx = 0.38381og[/^(2 keV)]/[/^ (2500 A)]. 



most data points at frequencies less than 10^^ Hz, in- 
cluding the upper limit in the radio. We integrate the 
radio-quiet template from 6 cm to 1 keV and the Chan- 
dra data from 1 keV to 10 keV and estimate the bolo- 
metric luminosity based on the full extent of our data to 
be Lboi(SED) w 1.3 x lO'*'^ ergs s~^, which is consistent 
with the value estimated from the UV correction. We 
take this value as our best estimate of Lboi- 

7.4. Black Hole Mass 

For Seyfert 1 galaxies and quasars, there exists a 
correlation between the radius of the broad-line region 
(^blr) and the luminosity of the optical non-stellar 
cont inuum (most often measured at 5100A), tblr, oc 
L^. iKaspi et al.l (|2000D derived this relationship based 
on their observations of 17 reverberation-mapped PG 
quasars combined with the results of 17 lower-luminosity 
AGNs with reverberat ion data. This r elationship was 
recently updated by Bentz et al.l (|2006D by taking into 
account the host galaxy starlight contamination. Us- 
ing this relation in tandem with the virial relation, 
Mbh = f 7'BLR A V^/G, one can estimate the mass of 
a black hole with a single spectrum. However, the nor- 
malization and slope in the radius-luminosity relation, 
and the value of the normalization factor / in the virial 
relation, which depends on the geometry and kinematics 
of the broad line region, are still subject to considerable 
uncertainty. In addition, different permutations of this 
relation not only use different emission lines for the line 
width measurement AV^, but they also measure AV^ in 
different way s, either using the F WHM or the line disper- 
sion (aiine)- lOnken et al.l (|2004l ) calibrated the average 
normalization factor, / = 1.4 by assuming all AGNs fol- 
low the same Mbh — f* relation as quiescent galaxies. 

We rely on the previous measurement of the H/3 
line width f rom a high-resolution Keck spectrum by 
iBarth et all (|2004D . They fit the broad H/3 line with 
a Gaussian model and found wfwhm = 765 km s~^. 
The conversion from FWHM to criino for a Gaussian fit 
is CTiinc ^ FWHM/2.35. 

Using t he radius-luminosity relation of 
IKaspi et al.l (|2000l ) and the broad H/3 li ne width of 
ypvyHM = 765 km s~^ discussed above, iBarth et al.l 
()2004l ) previously estimated the black hole mass of 
POX 52 to be Mbh « 1.6 x 10^ Mq. We start our 
new calculation of black hole mass using the combined 
revisions of Bentz et al. and Onken et al., and our new 
HST measurements, such that 



Mbh = 1.05x10^ 



A£a(5100 A) 
1044 ergs s-i 



0.518 



^FWHM 

10''' km s- 



(2) 

We note that our HST imaging is non-simultaneous with 
the previous Keck spectrum from Barth et al., and there- 
fore potentially introduces some error in the Mbh calcu- 
lation due to source variability. The extinction-corrected 
point source absolute magnitude calculated from the 
GALFIT decomposition is MF435vy(PSF) = -15.7. For 
comparison, we also calculate the extinction-corrected 
absolute magnitude needed to maximize the PSF con- 
tribution to the galaxy, Mi?435vi/(PSFMAx) = — 
16.3. These two magnitudes correspond to a range of 
ALa(5100 A) (3.2- 5.6) x lO^^ ergs s'^ Using wfwhm 



765 km s ^ from IBarth et al.l (|2004D . we calculate 
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TBLR = 2.0 - 2.7 It. days and A/bh » (3.1 - 4.2) x 10^ 

Mq. 

As noted earlier, there are many different recipes to 
estimate black hole mass that rely on different ver- 
sions of the r adius-luminosity and virial mass relations. 
iMcGill et al.l (j2008. ) recently compiled a list of some of 
the most common versions and we use their formula- 
tions to explore the range of possible Mbh values in 
POX 52. We eliminate those relations in which the slope 
of the radius-luminosity relation is > 5.5, which con - 
tradicts the updated measurement of lBentz et al.l ()2006D 
and are left with three alternate relations (equations 8, 
10 and 13 of McGill et al.). We use the same parame- 
ters above from o u r HS T measurements and the wfwhm 
from iBarth et al.l ()2004[ ) , and find a range of black hole 
masses of Mbh ~ (1-8 - 2.7) x 10^ M© for POX 52, 
with an average of (Mbh) = 2.2 x 10^ Mq. This then 
extends our total range of possible black hole mass to 
(1.8 - 4.2) X 10^ Mq, increasing the evidence that POX 
52 contains an intermediate-mass black hole. This mass 
range is due to the scatter in using slightly different rela- 
tions in the calculation of Mbh, but the sys tematic scat- 
ter in using any single relation is ^ 0.5 dex (jOnken et al.l 
120041: IVestergaard fc Petersonll2006l ). 

7.4.1. Black Hole Mass Estimates from X-ray 
Variability 

Both vb and excess variance have been shown to 
be ro ughly correlated with black hole mass (|Papadakisl 
|2004[ ). We investigate both of these methods in or- 
der to independently confirm our mass estimates from 
optical data. We first esti mate the expected z^g f rom 
the properties of the AGN. iMarkowitz et all (|2003f) fit 
a simple relation with I/vb measured in days, 1/vb = 
Mbh/10^'^ M0, in which vb is solely dependent on 
black hole mass. This relation seems to hold quite 
well for normal Seyfert 1 galaxies, but it underesti- 
mates black hole masses calculated from vb for NLSl 
galaxies, which have been known to have shorter break 
timescales (and therefore higher break frequencies) than 
classical Seyfert 1 galaxies with similar black hole masses 
(|Papadakisil2004D. In order to try to account for this 
disparitv. iMcHardv et al.l ()2006[ ) incorporated bolomet- 
ric luminosity in the vb-Mbh relation, in an effort to ac- 
count for the higher accretion rates seen in NLSl galaxies 
compared to the average Seyfert 1. Taking the inverse of 
Vb = ^/Tbi where Tb is known as the break timescale, 
McHardy et al. found the best fit in the Ts-MBH-^boi 
plane for their sample of Seyfert 1, NLSl, and Galactic 
black holes to be logTe = 2.101ogMBH - 0.981ogLboi - 
2.32, where Tb is measured in days, Mbh is measured 
in units of 10^ Mq, and Lboi is measured in units of 
IQii gj.gg g-i^ Q^j. average estimate of black 

hole mass and bolometric luminosity to calculate the 
corresponding vb from each of these relations and find 
i/B sa 1.2 X lO^"' Hz using the Markowitz et al. prescrip- 
tion and i^B ~ 3.8 x 10~^ Hz using the McHardy et al. 
prescription. Since our PSD analysis was unable to deter- 
mine Vb, we note that vb estimated from the Markowitz 
et al. relation occurs near the visual break seen in Figure 
[T3l and vb estimated from the McHardy et al. relation 
occurs at the high- frequency end of this same plot. 

Other recent work has explored alternative methods to 
correlate the variability in observations too short for PSD 



analysis with black hole mass. Nikolajuk et al. (2004, 
2006) calibrated the relationship between the excess vari- 
ance of the variability in a given observation and the mass 
of the black hole. The excess variance must be sampled 
from the high-frequency end of the PSD, above vb where 
the shape of the PSD is described by P{v) oc v'^. The 
cornerstone of this method is the assumption that the 
amplitude of the PSD ai vb multiplied hy vb is con- 
stant for all objects and independent of black hole mass, 
i.e. P{vb)i'b — constant. This has been shown to be 
roughly accurate for a handful of objects studied to date 
(IMarkowitz et al.ll2003HMcHardv etalll2004l : IPapadakid 
|2004[ ) . Since the excess variance is simply the integral of 
the PSD, it can then be related to vb and consequently 
AfBH- This technique has only been tested on a few low- 
mass black holes in other AGNs, namely NGC 4395 and 
the NLSl NGC 4051. 

Without a measurement oi vb for POX 52, we cannot 
determine that we are sampling the appropriate section 
of the PSD, and therefore, we continue on with the caveat 
that the mass estimates calculated using this method are 
tenuous at best. The Chandra observation had a higher 
signal-to-noise ratio than the XMM-Newton and there- 
fore had a better chance at producing a reasonable result. 
We are unable to probe timescales shorter than the break 
timescale predicted by the iMcHardy et al. (2006) rela- 
tion, Tb « 279 s, so we continue wi th this method using 
the simpler IMarkowitz et al.l (|2003f ) relation which gives 
a wide range of break timescales (Tb w 5460 — 13, 660 
s) corresponding to our full Mbh range. We tested the 
method using time bins of 100 s and 500 s and noted no 
significant difference between the two, except for the fact 
that using a smaller time bin allows us to divide the ob- 
servation into a larger number of segments and therefore 
probe lower break timescales. We calculated black hole 
masses using segment lengths of T « 2, 775 — 12, 500 
s with 100 s time bins and segments of length, Tb ~ 
6, 250 — 12, 500 s with 500 s time bins, resulting in an 
average black hole mass of A/bh = (3.9 ± 1.7) x 10^ Mq, 
where the errors represent the range of A/bh calculated 
for various Tb. If the break frequency of POX 52 is 
Vb ~ lO^'' Hz or lower, then we are able to probe the 
region above this frequency and calculate a black hole 
mass of Mbh ~ 3.9 x 10^ Mq. But if the break fre- 
quenc y is higher, near th e vb ~ 10~^ Hz predicted by 
the M cHardv et al.l ()2006D relation, we are unable to ex- 
plore the region of interest in frequency space and our 
black hole mass calculations are inconclusive. 

7.5. Accretion Luminosity 

Using the range of black hole masses we have estimated 
here, the Eddington luminosity of POX 52 is in the range 
L^dd = (2.9 — 5.5) X 10'''^ ergs s~^. Taking the bolometric 
luminosity calculated from the SED, we find an estimate 
of the Eddington ratio of /jboi/^Edd = 0.2 — 0.5. 

NLSl galaxies have been shown to have higher Edding- 
ton ratios on average than their Seyfert 1 counterparts, 
usually with L/L^dd ~ 1 for NLSl gala xies compared t o 
L/L^dd « 0.1 for a normal Seyfert 1 (jXu et al.l[2007| ). 
However, there are a few objects, e.g. NGC 4395, that 
strongly contradict this dichotomy. Despite nearly iden- 
tical spectra in the optical, POX 52 has an Eddington ra- 
tio that is more than two orders of magnit ude larger than 
the iboi/iEdd = 10"^ for NGC 4395 (jPeterson et all 
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l2005f) , demonstrating that there must be substantial dif- 
ferences between the central engines of these prototypical 
low-mass AGNs. An Eddington ratio of ^ 0.5 suggests 
POX 52 is not as highly accreting as many NLSl galaxies 
and may represent an intermediate object in the contin- 
uum between NLSl and Seyfert 1 galaxies. 

The unabsorbed spectral slope (F) of POX 52 remained 
relatively unchanged between X-ray observations, sug- 
gesting little ch a nge in the intrinsic continuum emission. 
IShemmer et al.l (|2006| ) studied the correlation between 
X-ray spectral slope and accretion rate (iboi/-^Edd) for 
AGNs up to z w 2, increasing the luminosity range to 
more than 3 orders of magnitude, almost twice as large 
as previous studies (e.g. Porquet et al. 2004, Wang et 
al. 2004, Bian 2005). Despite having a much lower mass 
than any of sources used in the Shemmer et al. study 
(Mbh« 10^ - 10^° Mq), pox 52 falls among the ma- 
jority of objects in the sample when plotting F versus 
Lhoi/ L-Edd-: making it consistent with this correlation, 
and suggests this relationship is independent of black 
hole mass. With F = 1.8 and iboi/-^Edd ~ 0.2 — 0.5, 
POX 52 falls in the overlapping region of the appar- 
ent distribution split between NLSl (F « 1.75 — 2.6, 
Lhoi/LEdd ~ 0.5 - 1.5) and Seyfert 1 (F w 1.5 - 2.3, 
iboi/^Edd ~ 0.05 — 0.5) galaxies. It is important to note 
that bolometric luminosity is most often estimated from 
a single-band using a bolometric correction (as done in 
Shemmer et al.), and that there might be significantly 
more scatter or bias in these relations because of these 
corrections. However, since Lboi for POX 52 is calcu- 
lated using a compilation of datasets over a wide range 
of wavelengths, it lends itself to a more secure calculation 
of the accretion luminosity of this object. 

7.6. Black Hole - Host Galaxy Correlations 

Relationships connecting the black hole mass with 
host galaxy properties, such as velocity dispersion, bulge 
mass, and bulge luminosity can be useful in predicting 
the mass of a black hole of many objects from single pho- 
tometric or spectroscopic observations. Determining the 
slope and amount of scatter in these correlations, specif- 
ically in the low-mass and high-mass ends of the black 
hole mass spectrum, can help us learn how black holes 
and host galaxies, evolve together in time (e.g. di Matteo 
et al. 2005). 

If we extra p olate the Mbh — cr* relation of 
iTremaine et al.l (l2002f) down to ct, = 36 km s^^ 
(jBarth et all \2004) . we find a black hole mass of 
Mbh = (1-4 ± 1.1) x 10^ M©, which is consistent with 
the average black hole mass calculated from the mass es- 
timators described above. We use the F435W magnitude 
of the outermost Sersic component in order to estimate 
an absolute i?-band magnitude for the host galaxy o f 
Mb « -16.0. Extrapolating the iMarconi fc HuntI (|2003D 
-^^BH~ibui relation down to a lower luminosity, we find a 
predicted Mbh of AfBH~ 8.7 x 10^ Mq, which is a factor 
of 2 — 5 hi gher than current estimates. The upd ated ver- 
sion o f the lMarconi fc HuntI (|2003,) relation from lGrahamI 
(|2007D predicts a black hole mass of Mbh~ 7.4 x 10^ Mq, 
which is a factor of 20 — 40 larger than estimated here. 

The measured host galaxy color places a con- 
straint on the age (and thereby mass-to-light ratio) 
of the galaxy stellar po pulations. The code kcorrect 
(|Blanton fc Roweiil2007f) fits a combination of templates 



with a range of star formation histories, metallicities, and 
reddening values and returns the stellar mass [assuming 
a Chabrier (2003) initial mass function]. Using either 
the HST measurements presented here, or the BVRI 
measurements from Barth ct al. (2004), we find a stellar 
mass of Mhost ^ 1-2 x 10^ Mq. If we extrapolate the 
iHaring fc Rb3 ()2004D Mbh - Mbui relation to the host 
galaxy mass of POX 52, the bulge mass implies a black 
hole mass of Mbh ~ 1.1 x 10^ Mq, a factor of 2 — 6 larger 
than our current estimate. 

The compact systems investigated in iGreene et al.l 
(l2008h follow the same trends. For the Mbh — ibui re- 
lation, the black hole masses are on average two orders 
of magnitude less massive than their host galaxy lumi- 
nosities would predict and are similarly an order of mag- 
nitude less massive than their host galaxy masses would 
predict. The divergent Mbh — -Mbui and Mbh — ^bui rela- 
tions in these galaxies as compared to elliptical galaxies 
may be an outcome of their different formation histo- 
ries, as reflected in their different fundamental plane loca- 
tions, which still leaves the possibility that the Mbh — ct* 
relation represents a tighter empirical correlation in this 
low-mass regime. The Mbh — relation seems to more 
accurately reflect the pr operties of lower ma ss black holes 
and their host galaxies (jBarth et al.ll2005l ). whereas the 
Mbh — ibui and Mbh — Mbui relations both overestimate 
the black hole mass by factors of ~ 2 — 100. 

There has been recent work suggesting a correlation 
between black hole mass and Sersic index for ellipti- 
cal galaxies or the bulges of disk galaxies, for the pur- 
pose of predicting black hole masses in other objects 
(|Graham fc Driver! l2007f ). The measured Sersic index 
n 4 and mass Mbh « 3 x 10^ Mq make POX 52 a 
strong outlier in this relation, which previously only in- 
cluded objects with Mbh « 10*^ -10^ Mq. For n = 4, the 
mass predicted by this relation is Mbh ~ 10* Mq, ap- 
proximately 3 orders of magnitude larger than all other 
mass estimates of POX 52 suggest. The inclusion of more 
objects like POX 52 and those in the GH sample demon- 
strate that there cannot be a single, tight Mbh—"- cor- 
relation applicable across the entire mass spectrum of 
galaxies. 

Recent HST imaging surveys have found evidence for 
a transition in the population of central star-cluster nu- 
clei in early-type galaxies, such that galaxies more lu- 
minous than Mb ~ —20.5 mag, or more massive than 
a few xlO^" Mq, tend to lack resolved central nuclei, 
while lower-lumi nosity early-type galaxies are more likely 
to be nucleated (Wehner fc Harris! 120061 : !Ferrarese et all 
!2006tlC6te et al.ll2006D . These groups further found that 
the black holes in luminous ellipticals and the stellar nu- 
clei in lower-mass galaxies followed the same relation- 
ship between central object mass and host galaxy mass. 
They suggested that formation of a central massive ob- 
ject might be a universal property of early-type galaxies, 
but that the central object would be either a massive 
black hole or a central star cluster, depending on the 
host galaxy mass. It has already been demonstrated that 
NGC 4395 has both a n AGN and a nuclear star cluster 
(|Matthews et al.l!l999! ) and our GALFIT decompositions 
suggest the possibility that POX 52 might as well. Obser- 
vations such as these provide further motivation to test 
for the presence of central black holes in dwarf ellipticals 
through dynamical searches or AGN surveys, and POX 
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52 is a unique object in this context, given its low stel- 
lar luminosity {Mb — —16.0 mag) and the unambiguous 
presence of a central black hole. 

8. SUMMARY AND CONCLUSIONS 

New HST observations have confirmed POX 52 to be 
a dwarf elliptical galaxy, with a Sersic index of n w 4.3, 
re ~ 1.15, and Mj = —17.3. These parameters, along 
with the previously me asured stellar velocity dispersion 
from' Barth et all ()2004f ). show POX 52 to reside near the 
dwarf elliptical family in the fundamental plane. POX 
52 was the first AGN found to reside in a dwarf elliptical 
galaxy, but more recent studies of low-mass AGNs and 
their host galaxies h ave found a number of similar objects 
(|Greene et al.l[2008h . 

We see strong variability over 500 — lO'' s timescales 
with both Chandra and XMM-Newton and a substantial 
change in the absorbing column density due to partial 
covering over a 9 month time period. Data from the 
Optical Monitor on XMM-Newton showed no significant 
variability in the UV over short ^ lO'^ — 10^ s timescales. 
We see no variability within the 2a errors in the 7 month 
period between the OM and GALEX NUV observations, 
but this does allow variations of 10 — 35% on the la 
uncertainty level. The Spitzer spectrum shows evidence 
of possible star-formation due to the PAH features and 
red continuum, but that analysis is deferred to a later 
paper. 

Including data from the VLA, Spitzer, GALEX, and 
2MASS, we compile an extensive SED that shows POX 
52 to be consistent with the radio-quiet template of more 
massive quasars and estimate a bolometric luminosity of 
Lhoi = 1-3 X lO^'^ ergs s~^ by incorporating all of the 
available data. Using a variety of mass estimators, we 
find a range of black hole mass Mbh ~ (1-8 — 4.2) x 10^ 
Mq. We find that we are unable to determine the break 
frequency from the Ghandra PSD, either due to a low 
signal-to-noise ratio or insufficient temporal sampling. If 
vb ~ 10^^ Hz, then our mass estimate from the X-ray 
excess variance of Mbh = (3.9 ± 1.7) x 10^ Mq is con- 
sistent with the optical mass estimates. Otherwise, new 
X-ray observations of POX 52 in an unabsorbed state 
with the larger effective area of XMM-Newton would be 
needed in order to determine vb and estimate the black 
hole mass with this method. 

We have shown POX 52 to have an SED consis- 
tent with a scaled-down version of more massive type 1 
quasars and that it is in a different accretion state than 



NGC 4395, as seen in the two orders of magnitude dispar- 
ity between Eddington ratios. The number of low-mass 
AGNs known is incre asing with new resul t s from large 
AGN surveys, such as iGreene "fcHol (I200I l2007bD . but 
POX 52 remains unique in that it is one of the nearest 
and best studied of the low-mass AGNs. 
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